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One of the most interesting and important phenomena predicted to occur in heavy-ion collisions is the local strong 
parity violation. In non-central collisions, it is expected to result in charge separation of produced particles along the 
system's orbital momentum. I report on results of the charge separation measurement in Au+Au and Cu+Cu collisions 
at v/sjviv — 200 and 62 GeV with the STAR detector at RHIC based on three-particle mixed harmonic azimuthal 
correlations. Systematic study of parity conserving (background) effects with existing heavy-ion event generators, and 
their possible contributions to the observed correlations are also present. 



1. Introduction 

A heavy-ion collision provides a unique environment 
to study particle interactions at a very high temperature 
and extreme density. Among the most important features 
of the system created in non-central heavy-ion collisions 
are the strong magnetic field B ~ 10 15 T (eB ~ 10 4 
MeV 2 ) [TJ [2j [3], and large orbital angular momentum 
L ~ 10 5 [4, 5 . Together with non-uniform particle den- 
sity and pressure gradient in the overlap area formed 
by the colliding nuclei (depicted by the color region in 
Fig. [T]), such initial conditions result in a variety of in- 
teresting physics phenomena which are currently under 
study at the Relativistic Heavy Ion Collider (RHIC). 
Among the most extensively studied is the anisotropic 

i i B - magnetic field 
L - orbital momentum 




reaction 
plane 



Figure 1: Schematic view of non-central heavy-ion collision. 
Colliding nuclei (depicted as dashed circles) are moving out- 
of-list. Magnetic field (B) and system orbital momentum (L) 
are perpendicular to the reaction plane (plane spanned by the 
impact parameter, b, and colliding nuclei direction). 

transverse flow, which reflects the evolution of initial 
space anisotropy of particle production in the overlap 
area into momentum space. This collective effect has 
been measured at RHIC and SPS (Super Proton Syn- 
chrotron) (see [6] and references therein), and the strong 
elliptic and directed flow observed at RHIC are key mea- 
sures which indicate the quark-gluon plasma formation 
in heavy- ion collisions. 



Another interesting phenomena predicted to occur in 
non-central heavy-ion collisions is global polarization and 
spin alignment [HE]. It manifests itself in preferential 
orientation of the spin of produced particles wrt. the 
system orbital momentum (perpendicular to the reaction 
plane as illustrated in Fig. [I]). This effect has been re- 
cently measured by the STAR Collaboration for strange 
hyperons (A, A) and vector mesons (if* , </>), and is ex- 
perimentally found to be consistent with zero within ex- 
perimental uncertainties [HJ [9] . 

In this presentation I report on results of the STAR 
studies of the local strong parity (V) violation. Experi- 
mental evidence for this effect would be an observation 
of charge separation along the direction of the magnetic 
field (preferential emission of the same charged particles 
in the same direction). This effect is illustrated by "+" 
(positive) and "-" (negative) charge separation in Fig. [I] 

The underlying physics of charge separation originates 
in the fundamental features of the QCD vacuum. The 
modern understanding of the QCD theory is that its 
vacuum (gluonic field energy) is periodic vs. so-called 
Chern-Simons number Nqs (related to the topological 
charge, Q, as Q = N C s(+oo) - Ncs(-oo)) [10 , and 
there exist localized in space and time solutions (topo- 
logical configurations) which corresponds to the trans- 
formation from one local minima to another. Transition 
between different vacua can be either via tunneling (in- 
stanton) or go-over-barrier (sphaleron) (see [3 and ref- 
erences therein). Figure [2] schematically illustrates these 
transitions. Quark interaction with either of these topo- 
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Figure 2: Schematic view of QCD vacuum and its topological 
transitions from one local minima to another (tunneling via 
instanton, or go-over-barrier via sphaleron). 
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logical configurations changes chirality, which is a V and 
time reversal odd transition. 

At the same time it is experimentally known that 
V and CP invar iance are (globally) preserved in strong 
interactions. The evidence comes form neutron electric 
dipole moment experiments, which put a strong con- 
strains on the possible magnitude of the parity violat- 
ing parameter [Hj [12]: < 1(T 10 (if ^ 0, then QCD 
vacuum breaks V and CP symmetry). 

It was recently argued by Kharzeev et. al. [U HJ |3] 
that in non-central heavy ion collisions (local) formation 
of metastable V-odd domains is not forbidden. Such lo- 
cal strong parity violation can be viewed as a combined 
result of the following physics effects [3 : 

• External magnetic field aligns quark spins along or 
opposite to its direction. Right-handed quark mo- 
mentum becomes opposite to the left-handed one. 

• Vacuum topological transitions produce a domain 
with local excess of left or right handed quarks 
Meft 7^ bright (what corresponds to non-zero topo- 
logical charge). 

• Magnetic field induces electric field which is parallel 
to it: E ~ ■ B. In this electric field positive 
and negative charges start to move opposite to each 
other, which in a finite volume results in charge 
separation. 

Depending on the sign of the domain's topological charge, 
positively charged particles will be preferentially emitted 
either along, or in the direction opposite to, the system 
orbital angular momentum, with negative particles flow- 
ing oppositely to the positive particles. 



2. Experimental observable 

Charge separation can be described by adding V-odd 
sine terms to the Fourier decomposition of the particle 
azimuthal distribution with respect to the reaction plane 
angle, 9rp [13]: 

dN °° 

— — oc 1 + 2 y v n ± cos nA(j) a + 2a± sin A(j) a , (1) 

' /0 - , . ' 

where A<p a = <p a — ^rp is the particle azimuth relative to 
the reaction plane, v n are coefficients accounting for the 
anisotropic flow (vi is called directed, and V2 is elliptic 
flow). The a parameters, a_ = — a+, describe the V- 
violating effects (in the present study we consider only 
the first harmonic). 

For mid-central collisions at RHIC energies, predicted 
value for the a parameter is of the order of 1% [1 , which 
is too small to be observed in a single event. The direct, 
V-odd, observable (a±) (average asymmetry over many 
events) must also yield zero, and the only way to exper- 
imentally probe charge asymmetry is to use multi- (two 



or more) particle correlations, what limits the measure- 
ment to V-even observables only. This brings a problem 
of separating (or suppressing) ^-conserving background 
physics effects in the measurement. Physics backgrounds 
can be sorted into two different categories: background 
which produce correlations independent of the reaction 
plane orientation and those which are reaction plane de- 
pendent. To suppress reaction plane independent back- 
grounds, Voloshin proposed to use the following two par- 
ticle correlator wrt. the reaction plane [T3] : 

(cos(0 a + 0(3 - 2V RP )) = (2) 
= (cos Acj) a cos A<j)p) — (sin Acj) a sin Ac/) p) (3) 

= [(vi^Vl^) + B in] - [(daap) + B out ]. (4) 

Here Bi n and B out correspond to the in- and out-of plane 
background contributions to the correlator, and a, (3 de- 
note particle charge. Figure [3] illustrates the relation be- 
tween various azimuthal angles in the laboratory frame 
and the reaction plane orientation. 




Figure 3: Relation between azimuthal angles of the particle 
a and (3 in the laboratory frame and the reaction plane. 

The correlator (|3| represents the difference between 
correlations projected onto an axis in the reaction plane 
and the correlations projected onto an axis that is out- 
of-plane or perpendicular to the reaction plane. The key 
advantage of using a difference equation is that it removes 
all the background correlations among particles a and f3 
that are not related to the reaction plane orientation [14], 

US]. 

3. Expectations for charge correlations 

The first estimates [1 of the charge asymmetry pre- 
dicted a signal of the order of \a\ ~ Q/N n +, where 
Q = 0, =bl, ±2, ... is the topological charge and N n + is 
the positive pion multiplicity in one unit of rapidity. 
More accurate estimates [2 were found to be close to 
the same number, of the order of 10 -2 for mid-central 
collisions, which corresponds to 10 -4 for the correlator 
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(aaap) (third term in Eq.j4|. Without in medium effects, 
one expects the symmetry between same- and opposite- 
sign correlations: (a + a + ) = (a_a_) = — (a + a_) > 0, 
but in case of heavy ion collisions one needs to account 
for correlation modification due to particle interaction 
with the medium. This is similar to the back-to-back 
suppression of jet-like correlations, and is predicted to 
result in in suppression of opposite sign correlations [2]: 
(a + a + ) = (a_a_) ^> — (a+a_). The correlator (a a ap) 
should follow (a typical for any kind of correlations due 
to clusters) 1/N c h dependence (N c h is the charge mul- 
tiplicity). Local parity violation is a non-perturbative 
effect and the main contribution should come from par- 
ticles which have transverse momentum smaller than 
1 GeV/c [2 , with correlation width of the order of one 
unit in rapidity (typical size of the V-odd domain is about 
1 fm). 

4. Measurement technique 

The reaction plane orientation is not known a priori 
in the experiment, and it needs to be reconstructed from 
particle azimuthal distributions [16]. For that reason in- 
stead of measuring the two particle correlator wrt. the 
reaction plane ([2]), we evaluate the three particle correla- 
tor, where the third particle, c, gives an estimate of the 
reaction plane. Such a three particle correlator can be 
related to the correlator ([2| by the following equation: 

(cos(0 a + (j)? - 2& RP )) = (cos(0 a + fa - 2<j> c ))/v 2 , c - (5) 

Here V2 jC is an elliptic flow value of the third particle c, 
and it accounts for the finite resolution in our estimate of 
the reaction plane angle [15j [16j [17] . Equation [5] assumes 
that particle c is correlated with particles a and (3 only 
via common correlation to the reaction plane. This third 
particle further complicates experimental data analysis, 
since it potentially adds to the measured values of the 
correlator ([5| some background contributions from gen- 
uine 3-particle correlations. These contributions are dis- 
cussed in Sec. El 

5. Experimental setup 

The results presented here are based on 14. 7M Au+Au 
and 13. 9M Cu+Cu collisions at the incident energies 
y/s]^=200 GeV, and 2.4M Au+Au and 6.3M Cu+Cu 
events at ^/s~nn=62 GeV. Collisions were recorded with 
the STAR detector during the 2004 and 2005 runs. 
Charged particle tracks were reconstructed in a cylindri- 
cal and azimuthally symmetric Time Projection Cham- 
ber (TPC) [18j [19]. A minimum bias trigger was used 
during data-taking with offline collision vertex cut of 
30 cm along the beam line from the center of the main 
TPC. 

The correlations are reported in the pseudorapidity re- 
gion 1 77 1 < 1.0 with particle momentum 0.15 < pt < 



2.0 GeV/c (unless stated otherwise). Standard STAR 
track quality cuts are applied [20] . The Large elliptic flow 
measured at RHIC [21] is used to estimate the reaction 
plane from particle distributions in the main TPC and 
two Forward Time Projection Chambers (FTPC) [22] . 
The latter cover pseudorapidity intervals 2.7 < \rj\ < 3.9. 
In the most forward direction, STAR has two Zero De- 
gree Calorimeter - Shower Maximum Detectors (ZDC- 
SMD) [231 [24]. ZDC-SMDs are sensitive to the directed 
flow of neutrons in the beam rapidity regions, which is 
used for the reaction plane reconstruction. 



6. Understanding detector effects 

Detector effects have been corrected with recentering 
procedure [16]. Corrections are applied run- by-run, sep- 
arately for positive and negative particles, for each cen- 
trality bin, and as a function of particle pseudorapidity 
and transverse momentum. The validity of the recen- 
tering method was verified by calculating three-particle 
cumulants according to [T4j [25] . 

Figure [4] shows three-particle correlator (cos(</> a + 0/3 — 
2(j) c )) as a function of reference multiplicity (measured 
charge multiplicity in \rj\ < 1) in Au+Au collisions at 
\/snn=200 GeV (a) before and (b) after correcting for 
detector effects. The signal is scaled by the reference 
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Figure 4: (cos(<p a + 4>p — 2<p c )) as a function of reference 
multiplicity for different charge combinations [20]: (a) be- 
fore corrections for acceptance effects, and (b) after correc- 
tions. In the legend the signs indicate the charge of particles 
a, /3, and c. 

multiplicity to better indicate the acceptance effects for 
central collisions. All the differences which do not de- 
pend on the relative charge of a and (3 disappear after 
the acceptance correction. We have performed several 
additional checks to ensure that the signal is not due to 
detector effects: 

• Distortions in the track momenta due to the charge 
buildup in the STAR Time Projection Chamber 
(TPC) at high accelerator luminosity. Results ob- 
tained from RHIC Run II (a low luminosity run), 
and those from Run IV divided into high and low 
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luminosity events, yield the same signal within sta- 
tistical uncertainties. 

• Signal dependence on reconstructed position of 
the collision vertex along the beam line has been 
checked, and no dependence has been found. 

• Displacement of track hits when it passes the TPC 
central membrane. The correlator ([5| has been cal- 
culated using only particles with the entire track in 
one of the different half-barrels of the TPC. Cor- 
rected for the signal dependence on the track sep- 
aration in pseudorapidity, results are found to be 
consistent with those obtained before introducing 
the rapidity separation. 

• Feed-down effects from non-primary tracks (i.e. 
resonance decay daughters) have been studied via 
cuts on track distance of closet approach (dca). Re- 
sults for dca < 1 cm and dca < 3 cm are found to 
be consistent within statistical errors. 

• Electron contribution to the measured signal has 
been checked via specific energy loss (dE/dx) in 
the volume of the TPC and found to be negligible. 
This check was done to test whether the signal was 
due to hadron production, or lepton production. 

• Studied a correlator similar to ([5| but with the re- 
action plane angle rotated by 7r/4. This new corre- 
lator should only deviate from zero due to detector 
effects. It was found to be consistent with zero 
within statistical errors. 

• Variation depending on the charge of the third par- 
ticle used to reconstruct the reaction plane and 
changes of the STAR magnetic field polarity. The 
variations does not change the observed signal. 

Our conclusion from a number of tests performed is 
that detector effects are not responsible for the observed 
correlations. 



7. Testing sensitivity to 2-particle 
correlations wrt. the reaction plane 

Figure [5|a) compares the three-particle correlations 
(cos((/) a + — 2(j) c )) obtained for different charge com- 
binations, as a function of centrality, when the third 
particle is selected from the main TPC (solid sym- 
bols) with when it is selected from the Forward TPCs 
(open symbols). Elliptic flow data have been taken from 

Ref. [251 E2- 

Assuming that the second harmonic of the third par- 
ticle is correlated with the first harmonic of the first two 
particles via a common correlation to the reaction plane, 
the correlator should then be proportional to the elliptic 
flow of the third particle. Figure [5jb) shows very good 
agreement between the same charge correlations obtained 
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Figure 5: (a) (cos((j> a + cjyp — 2(j> c )) obtained by selecting the 
third particle from the main TPC and the Forward TPCs [20] . 
(b) (cos(4> a + 0/3 — 20 c )) after scaling by the flow of the third 
particle. The shaded areas represent the systematic uncer- 
tainty from V2, c scaling. 



with the third particle in the TPC and FTPC regions, 
which supports for such correlations the assumption of 
Eq. [5] The opposite charge correlations are small in 
magnitude and with current statistics it is difficult to 
conclude on validity of the same assumption. Results 
obtained with the event plane reconstructed with ZDC- 
SMD are consistent with those shown in Fig.[5|b), though 
the statistical errors on ZDC-SMD results are about 5 
times larger. 

Agreement between TPC, FTPC, and ZDC SMD re- 
sults give us a confidence that 3-particle correlation that 
we measure are indeed sensitive to the 2-particle correla- 
tions wrt. the reaction plane. 



8. Modeling physics backgrounds 

Figure [6] shows the correlator ([5| for Au+Au colli- 
sions at v AsjVAr=200 GeV. Blue symbols and lines mark 
opposite-charge correlations, and red are same-charge. 
Star symbols represent the STAR data, with shaded area 
reflecting the uncertainty due to V2 scaling. Other sym- 
bols in Fig. [6] show possible background two particle cor- 
relations wrt. the reaction plane estimated from various 
Monte- Carlo event generator with true reaction plane: 

• HIJING [28] with default, quenching-off setting: 
with (triangles), and without (diamonds) V2- 
afterburner. Elliptic flow is added by "shifting" 
method [16 , with V2 values consistent with STAR 
measurements at the given centrality. 

• UrQMD [29] with default setting (circles). 

• MEVSIM [30] Monte-Carlo simulations (squares). 
MEVSIM tests resonance (</>, Z\, p, cj, and K*) con- 
tribution with realistic elliptic flow pattern. 

Thick solid lines in Fig. [6] indicate genuine three particle 
correlations from HIJING (corresponding estimates from 
UrQMD are about factor of two smaller). HIJING three 
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Figure 6: Comparing correlator ^ from STAR data [20] 
vs. various physics background event generators in Au+Au 
collisions at ^snn = 200 GeV. None of the generators can 
reproduce same-charge correlations. 



particle correlations produce data-like opposite-sign sig- 
nal (compare blue line with blue stars in Fig. [6|, which 
may indicate dilution of the opposite-sign signal with ef- 
fects not related to the reaction plane orientation. 

Figure [6] shows that no generator gives qualitative 
agreement with the data, thought they produce non-zero 
correlations. These models do not match the correlations 
for (cos((j) a — 4>p)) that are seen in the data either, which 
points to the need for better modeling of two-particle cor- 
relations to give quantitatively meaningful comparisons 
for correlator ([5|. Other background effects studied are 
(see for details^Q]): 

• Correlations from processes in which particles a 
and /3 are products of a cluster (e.g. resonance, 
jet, di-jets) decay, and the cluster itself exhibits el- 
liptic flow or decays (fragments) differently when 
emitted in-plane compared to out-of-plane 

• Jets as potential source of reaction-plane dependent 
background since their properties may vary with 
respect to the reaction plane. 

• Cluster formation, which plays an important role in 
multiparticle production at high energies, and may 
account for production of a significant fraction of 
all particles. 

• Directed flow, which on average is zero in a sym- 
metric pseudorapidity interval, but can contribute 
to the correlator J5| via flow fluctuations (see the 
first term in Eq. [if7 

• Effect of global polarization (discussed in the intro- 
duction), which produce polarized secondary par- 



ticles along the direction of the system's angular 
momentum. 

From our studies we conclude that none of the back- 
ground effects listed above can be responsible for the 
same-sign correlations shown in Fig. [6] 



9. Results 

Figure [7] shows correlator ([5| in Au+Au and Cu+Cu 
collisions at ^snn — 200 GeV. The correlations are 
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Figure 7: Correlator ([51 in Au+Au and Cu+Cu collisions at 
■^/snn — 200 GeV [20f . The shaded band indicates uncer- 
tainty associated with v 2 scaling. Thick solid (Au+Au) and 
dashed (Cu+Cu) lines represent HIJING estimates for pos- 
sible non-reaction-plane dependent contribution from many- 
particle correlations. 

weaker in more central collisions compared to more pe- 
ripheral collisions, which partially can be attributed to 
dilution of correlations which occurs in the case of par- 
ticle production from multiple sources. The correlations 
in Cu+Cu collisions (open symbols), appear to be larger 
than the correlations in Au+Au (solid symbols) for the 
same centrality of the collision. One reason for this dif- 
ference may be the difference in number of participants 
(or charge multiplicity) in Au+Au and Cu+Cu collisions 
at the same centrality. The local ^-violation signal is 
expected to have approximately l/N c h dependence, and 
at the same centrality of the collision the multiplicity 
is smaller in Cu+Cu collisions compared to Au+Au. 
The difference in magnitude between same and oppo- 
site charge correlations is considerably smaller in Cu+Cu 
than in Au+Au, qualitatively in agreement with the sce- 
nario of stronger suppression of the back-to-back corre- 
lations in Au+Au collisions. 

Figure [8^ a) shows the dependence of the correlator ^ 
on the difference in pseudorapidities of two particles, 
\Vat ~ Vpli f° r 30-50% centralities, while Fig. [8^b) shows 
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Figure 8: Au+Au at 200 GeV [20]. The correlations depen- 
dence on (a) pseudorapidity separation, \rj a — ^7/3 1 , and (b) 
on (pt,a +pt,p)/2 for centrality 30-50%. The shaded band 
indicates uncertainty associated with V2 scaling. 



its dependence on the sum of the transverse momentum 
magnitudes (without upper p t cut) for the same central- 
ity region. The correlations have a typical hadronic width 
of about one unit of pseudorapidity, but we do not ob- 
serve them to be concentrated in the low p t region as 
naively might be expected for V- violation effects. 

Other results, in particular for collisions at 
y / S7V7v=62.4 GeV, correlation scaling with number 
of participants, and their dependence on transverse 
momentum difference, can be found in [20] . 



10. Summary 

Formation of local P-odd domains has been predicted 
in nuclear collisions [1 . This effect should result in 
charge separation along the orbital momentum of the 
system created in non-central heavy-ion collisions. The 
STAR measurements of 3-particle azimuthal correlations 
in Au+Au and Cu+Cu collisions at v AsjVAr=200 and 
62 GeV reveal non-zero signal. Measured 3-particle az- 
imuthal correlations are directly sensitive to local strong 
parity violation, but susceptible to contributions from 
^-conserving backgrounds. So far we could not explain 
the observed same sign correlations, and data can not 
be describe with any of the existing Monte-Carlo event 
generators of the heavy- ion collisions. At the same time, 
qualitatively the data agrees with predictions for local 
^-violation (though the signal persists to higher trans- 
verse momentum than expected). The presented mea- 
surements demand detailed theoretical calculations of the 
^-violating signal and backgrounds to make a defini- 
tive statement about possible local strong P-violation in 
heavy ion collisions. 
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